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ABSTRACT 

We analyze the optical, UV, and X-ray microlensing variability of the lensed quasar SDSS 
J0924-I-0219 using six epochs of Chandra data in two energy bands (spanning 0.4-8.0 keV, 
or 1-20 keV in the quasar rest frame), 10 epochs of F275W (rest-frame 1089A) Hubble Space 
Telescope data, and high-cadence it-band (rest-frame 2770A) monitoring spanning eleven years. 
Our joint analysis provides robust constraints on the extent of the X-ray continuum emission 
region and the projected area of the accretion disk. The best-fit half-light radius of the soft 
X-ray continuum emission region is between 5 x 10^^ and 10^^ cm, and we find an upper limit 
of 10^® cm for the hard X-rays. The best-fit soft-band size is about 13 times smaller than the 
optical size, and roughly TGMbh/ c^ for a 2.8 x IO^Mq black hole, similar to the results for other 
systems. We find that the UV emitting region falls in between the optical and X-ray emitting 
regions at 10^'^ cm < ri/ 2 ,uv < 3 x 10^® cm. Finally, the optical size is significantly larger, 
by l.Scr, than the theoretical thin-disk estimate based on the observed, magnification-corrected 
/-band flux, suggesting a shallower temperature profile than expected for a standard disk. 


1. Introduction 


Quasars are now known to play a major role in 
galaxy formation and cosmology through their 
effects on star for mation and ionization (e.g., 
iTevssier et al. 2011 ). However, the detailed struc¬ 
ture of the central engine has remained unclear, 
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as well as the relationships between emission 
at different wavelengths. The main difficulty is 
that the luminous source is so compact, requiring 
sub-microarcsecond resolution for direct imaging. 
Studies have made use of the time domain to 
probe their inner structure, such as placing upper 
limits on the variable source siz e using character¬ 


istic t imescale arguments ( e.g.. ISmith fc Hoffleit 


19631 : I Vaughan et al.l 120031 ) or using reverbera¬ 


tion mapping to probe the broad line region and 
calib rate black hole masses (e.g., IPeterson et al 


19921) . Reverberation lags in highly variable AGN 
have been used to exp lore the continuum emission 
( Shappee et alj_2014ri and the X -ray generation 


process (e.e.. ICackett et ah 20131) . Variations ir 
the profiles of Fe fluorescence lines can also be 
used to constrain the inner and outer radii of the 
Fe lin e reflection component of the a ccretion disk 
le.g.. iFabian et allll989l : Laor 199lh . Microlens¬ 
ing of gravitational lensed quasars provide another 
way to “zoom in” on the most compact regions in 
AGN and dissect their structures. 
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Microlensing in strongly lensed quasars can 
be used to measure many physical properties 
of the quasar continuum source and the lens 
galaxy. Si nce quasar micro l ensing was first de- 
scr ibed by Chang fc Refsdall (1197911 and detected 
by ICorrigan et al.l ( 199lll . studies have success¬ 
fully isolated its signal by analyzin g observations 
of le n sed quasars over time (e.g., Morgan et al. 
2011)1 : Mosouera et al. 2013h. wav elength (e.g., 


2013t iRathna Kumar et al.l l2013f) . In particular. 


Bate et al .II2OO8L [Munoz et al.ll201l[l. or both (e .g., 


Poindexter et al. 20081 : Eigenbrod et al.ll2008J bl. 

Microlensing by stars in the lens galaxy can 
(de)magnify individual images in a multiply- 
imaged system, causing the ratios between the 
image fluxes to be significantly different from the 
predictions of strong lens models fit to the ob¬ 
served macroscopic properties of the lens. Due to 
relative motions between the source, lens, stars 
in the lens, and observer, these magnifications 
change with time, leading to uncorrelated vari¬ 
ability in the light curves of lensed quasars. 

The analysis of quasar flux ratio anomalies 
in single- or few-epoch multi-wavelength obser¬ 
vations is a relatively inexpensive means of es- 
timating the quasar continuum source size (e.g., 
Poolev et al. 2007t Bate et al. 2008t Flovd et al 


2009 : Mediavilla et al.ll2011 : Mosauera et al.l2()Tl ). 

While insensitive t o the mass function (e.g., 
Lewis fc Irwin 19951) . this requires assuming the 
mass of an average star in the lens galaxy (M) 
( Wambsganss et ahl Il992l) and that the macro¬ 
model lens magnifications and differential extinc¬ 
tions are well understood. The choice of a lin¬ 
ear or logarithmic prior on the source size also 
tends to significantly affect th e resulting size es¬ 
timates ( Mosauera et al. 201 ih . and the estimates 
can be further complicated for lenses with sig¬ 
nificant time delays, where the flux ratios can 
be contaminated by the intrinsic variability of 
the quasar. Moreover, microlensing magnifi¬ 
cation patterns are complex and frequently di¬ 
vided into “active” regions of many overlapping 
caustics and “passive” valleys with few caustics, 
and single-epoch data provide little constraint on 
the location of the source relative to these pat¬ 
terns. Mo nitoring lenses over time addresses these 
problems ( KochanekI 2004 : iKochanek et al.lbood 
Goicoechea et al.1 20081) while simultaneously per¬ 


mitting time de lay measurements (e.g.. recent 


the B ayesian Monte Carlo method of IKochanek 
( 20041 ) enables the extraction of probability den¬ 
sities for the physical parameters of interest (e.g., 
source size, (M), etc.) by modeling the microlens¬ 
ing variability using a Monte Carlo simulation of 
the observed light curves. 

A common conclusion of essentially all mi¬ 
crolensing studies is that there is a systematic dis¬ 
crepancy between microlensin g size measurements 
and t heoretical thin disk sizes (IShakura fc Sunvaev 
197.81 ) or the size estimated from the quasar’s 


ODse 

Blac 

rvea lummosiw le.g 
cbnrne et al. boil: 

. rooiev et ai.i zuunc 

VIosQuera & Kochanek 

2011 

: Munoz et al.1 2011: 

Hainline et al. 20131). 


The larger disk sizes inferred from microlens¬ 
ing suggest a shallower temperature profile for 
the disk, as might be found if irradiation of 
the outer disk by the i nner disk is important 
(see Morgan et al.1 20101) . It does not appear 


to be explainable by emission line contamination 
or sc attering (see Dai et al. 20ld Morgan et al 


2 OIOI ). The temperature profile can be deter¬ 


mined by measuring the dependence of the disk 


size on wavelength (e 

g.. Poindexter et al. 

2008; 

Eigenbrod et al.1 2008. 

al Blackburne et al.1 

2014; 

Jimenez-Vicente et al. 

2014). but there is a 

large 


scatter in the best-fit dependence among studies. 
The X-ray emission regions are muc h more 


2001 

Chartas et al. 

2006 

Poolev et al.1 



i et al.l2003: Blackburne et al 


ray source must be more compact. In fact, the 
analyses of microlens i ng var i ability for six quasars 


(Morgan et al. 

2008. 

2012: Chartas et al.1 

2012 

Mosauera et al. 

2013 

: Blackburne et al.l 

2013 


20141) from our Chandra monitoring program have 
consistently shown that the projection of the X- 
ray continuum emission region is confined to a 
region close to the inner edge of the accretion 
disk. The energy structure of the X-ray con¬ 
tinuum source can be explored by considering 
the hard and soft X-ray bands separately, and 
these studies have presented mixed evidence for 
an en ergy-dependent X -ray size. For example, 
while Chen et al. ( 201ll) found a more compact 


hard X -ray source in Q2237-l-0305. IChartas et al 


measurements bv lTewes et ahllioi^ lEulaers et al 


(I 2 OI 2 I) found that the softest X-rays were mi- 
crolensed more strongly than the hardest X-rays 
in image D by a factor 1.3 in RXJ1131 —1231. 
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Meanwhile, iBlackburne et al.l (120141) found no dif¬ 
ference between the hard and soft X-ray sizes in 
HE 0435-1223. 

The quad ruply-imaKed lens ed quasar SDSS 
J0924-I-0219 ( Inada et al. 2003 1. at a source red- 
shift of 1.524, shows the most anomalous optical 
flux ratio of any four-image lens, and it has been 
unclear whethe r the anomaly is du e to microlens- 
ing by stars (Peeples et al.l 2004 Keeton et al 


20061 : Morgan et al .1 200^ Flovd et al.l jpooir or to 
perturbations to the lens potential by dark matter 
substructure ( EigenbrpdeLalJ 2QQ6; FaureelM 


2011 : Bate et al. 201l[ "ochanek fc Palall 20041) . 
A comparison of the broad emission line flux 
ratios to the continuum flux ratios suggests 
that both substructure and microlensing must 
significantly contribute to the anomalous A/D 
flux ratio ( Keeton et al. 20061) . Furthermore, 
variability in the X-ray continuum flux ratios 
( Poolev et al. 2012 : Chen et al. 2012) also indi¬ 
cates tha t microlensing must be present in this 
system. [Morgan et al. ( 20061) analyzed the mi¬ 
crolensing using optical monitoring data from 
the SMARTS telescope and found that the ac¬ 
cretion disk radius as measured from microlens¬ 
ing is larger than, but marginally consistent 
with, that inferred from thin disk theory and 
the magnifi ca tion-c orrected /-band flux. Later, 
Flovd et al.l (|2009l ) found the size-wavelength 
power-law index R\ oc to be less than ( < 1.34 
for SDSS J0924-I-0219 at 95% confidence, challeng¬ 
ing the standard thin disk model (where C, = 4/3), 
although such single-epoc h studies general l y yiel d 
prior-dependent results ( Mosauera et al.l 120111) . 
Assuming that only microlen sing is respons i ble for 
the observed flux anoma ly, iPeeoles et al.l (l2004h 
and [Morgan et al] ( 2006h predicted a substantial 
brightening of the faint saddle point image D on 
a timescale of roughly one decade. 

Here, we extend the microlensing analysis of 
Morgan et al.l (20^ using recent optical, ultra¬ 
violet (UV), and X-ray monitoring observations 
of SDSS J0924-I-0219. While image D has re¬ 
mained roughly an order of magnitude fainter than 
image A, allowing for the possibility of millilens- 
ing by dark matter substructure, our light curve- 
based microlensing analysis remains robust. In 
Section[2l we describe the /?-band monitoring data 
obtained by the SMARTS and Euler telescopes, 
the UV data from the Hubble Space Telescope 


(HST), and the X-ray data from Chandra. In 
Section [3l we describe our analysis technique and 
present our results from a joint Bayesian Monte 
Carlo microlensing analysis. In Section 31 we 
discuss our results. Where needed, we adopt a 
Dm = 0.3, Da = 0.7 flat cosmology. 


2. Observations 


2.1. Optical Monitoring 

We monitored SDSS J0924-I-0219 for ten sea¬ 
sons (2003-2014) in the R band using the SMARTS 
1.3 m te lescope with the AN DICAM optical/infrared 
camera ( DePov et al.ll200'3 'n. and for three seasons 
(2010-2013) using the 1.2 m Euler Swiss Telescope 
as a part of the COSMO GRAIL projedH. In the 
rest frame of the quasar a,t z = 1.524, the effective 
wavelength of the R band corresponds to 2770A. 
To obtain light curves of the quasar ir nages, we use 
the ph otometric pipeline described in Tewes et al.l 
(l2013l) . which is based o n the deconvolution algo¬ 


rithm of iMagain et al.l (Il998l) . We fit four point 


sources plus a Sersic prof ile to the lens com p onent s 
using astrometry from Eigenbrod et al.l ( 2006 ). 
Unfortunately, we were unable to reliably sepa¬ 
rate the flux of image D from that of image A, so 
we will fit the combined flux for images A and D 
with our Monte Carlo models. A sample SMARTS 
image and its deconvolution are shown in Figure[I] 
along with the HST //-band and UV images for 
reference. We adjust the SMARTS and Euler 
light curves for images A-l-D, B, and C to mini¬ 
mize their dispersion in the overlapping observing 
seasons (2010/2011 and 2013). This adjustment 
compensates for any small color terms introduced 
by differences in the filters and quantum efficiency 
curves of the CCDs, and it corrects for any PSF 
variability caused by minor differences in the flat¬ 
ness of each detector. Our final optical light curves 
are shown in Figure [2] and provided in Table [T] 
(each “season” corresponds to a grouping of data 
points). The average seeing was I'.'S, with a range 
of 0"9 to 2"0. The absolute photometric calibra¬ 
tion of these light curves is obtained using SDSS 
DR9 r-band photometry of field stars. We esti¬ 
mate the zero-point uncertainty to be 0.07 mag. 
This calibration uncertainty is not included in the 


http: //www.astronomy.ohio-state.edu/ANDICAM/ 
http: //www.cosmograil.org/ 
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Fig. 1.— SMARTS images for SDSS J0924+0219 (top 
two panels), showing the raw data (left) and the deconvo¬ 
lution (right) from May 1, 2004. Bottom-left: HST H-bsLud 
image. Bottom-right; stacked HST UV image (ten epochs 
over two years). 


errors bars in Figure [2] or Tabled! as the absolute 
calibration has no effect on our further analysis 
which only uses magnitude differences. 

2.2. UV Data 

SDSS J0924+0219 was observed at 10 epochs 
between 2009 August and 2011 August using the 
Wide Field Camera 3 (WFC3) on board the Hub¬ 
ble Space Telescope (HST). The observations used 
the F275W filter in the UVIS channel, provid¬ 
ing flux ratios at 2750A (1089A in the source rest 
frame). The details of th e reduction of thes e obser - 
vations are prese nted in iBlac kburne et al. (12014 1. 


We u se imfitfits (IMcLeod et al.Nl998t iLehar et al 


2000 !) to perform PSF photometry and assign error 
bars that correspond to the difference between the 
flux ratio estimates from PSF and aperture pho¬ 
tometry. We then use the PHOTFLAM, PHOT- 
PLAM, and PHOTZPT header keywords to con¬ 
vert the instrumental magnitudes to the AB sys¬ 
tem (in this filter, the offset from ST to AB mag¬ 
nitudes is mAB — msT = 1.532). We provide the 
calibrated photometry in Tabled) We are unable 
to provide photometry for nearby comparison ob¬ 
jects due to the lack of other UV-bright sources in 
the small HST field. Figure [3| shows the final UV 


light curves. Note that while images A and D are 
merged in the optical data, they remain separate 
in the UV light curve analysis, with a mean flux 
ratio fA/fo = 20 ± 1. 

As the lens galaxy is expected to be faint in 
the UV, we searched the stacked UV image for the 
fifth, central image of the quasar which is expected 
from the lens model. To test for a central image, 
we used PSF fitting to carefully subtract the wings 
of the four bright quasar images. The PSF model 
was constructed from the exposure-weighted sum 
of ov ersampled and ap propriately rotated Tiny 


Tim (Krist et al.M2011h model PSFs for each of 


the 10 epochs. We fixed the relative positions of 
the five quasar images using the positions of the 
images and the lens galaxy from the CfA-Arizona 
Space Telescope Lens Survey (CASTLE^, allow¬ 
ing the overall position offset and the five normal¬ 
izations to vary. The residuals suggest the pres¬ 
ence of a fifth image, but the detectability is un¬ 
stable under changes in the PSF model. We can 
set an upper limit on its flux relative to image B of 
/s/Zb < 0.006. Such a detection, or limit, is not 
in conflict with the expectations from models of 
early-type galaxies (IKeeton et al.ll2005ll and could 
be further weakened by any dust at the core of the 
lens galaxy. Of course, a detection might also be 
problematic because it could also arise from nu¬ 
clear star formation or weak AGN activity in the 
lens galaxy, which would be difficult to distinguish 
from the source quasar emission. 

2.3. X-ray Data 

We monitored SDSS J0924-I-021 9 in X-rays us¬ 
ing th e ACIS imaging spectrometer (jCarmire et al. 
20031) on the Chandra X-Ray Observatory for six 
epochs between 2005 February 24 and 2010 Oc¬ 
tober 5. These observations were a component of 
a larger Chandra Cycle 11 monitori ng program. 


the de tails of which are published in IChen et al 
( 2012 !) . As with the UV data, images A and D 


remain separate in the X-ray light curve analy¬ 
sis. The full-band data are shown in Figure [3l 
Note that while we only show the relative X-ray 
m agnitudes, t h e cal ibrated fluxes are available 
in Chen et Tl ( 2012ll . We also use magnitudes 


rather than simply logarithms of the flux so that 
the X-ray flux ratios are on the same scale as the 


^http: //www.cfa.harvard.edu/castles/ 
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Table 1 

SDSS J0924+0219 Optical Light Curves 


HJD — 2,450,000 QSO A+D QSO B QSO C Source 


5543.806 

18.786 ±0.011 

20.243 ± 0.027 

21.720 ± 0.136 

SMARTS 

5544.828 

18.793 ± 0.003 

20.255 ± 0.007 

21.434 ± 0.018 

Euler 

5548.821 

18.794 ± 0.003 

20.258 ±0.007 

21.495 ±0.045 

Euler 

5550.755 

18.763 ±0.022 

20.291 ±0.073 

21.382 ±0.197 

SMARTS 

5558.831 

18.799 ± 0.004 

20.271 ±0.012 

21.370 ±0.024 

Euler 

5559.694 

18.782 ± 0.016 

20.253 ±0.029 

21.425 ±0.143 

SMARTS 

5564.789 

18.765 ± 0.003 

20.242 ± 0.009 

21.451 ±0.054 

Euler 

5565.761 

18.773 ±0.022 

20.262 ±0.077 

21.446 ±0.130 

SMARTS 

5572.758 

18.723 ±0.010 

20.170 ±0.047 

21.506 ±0.291 

SMARTS 

(5573.771) 

18.711 ±0.005 

20.179 ±0.007 

21.350 ±0.018 

Euler 



Note.— Light curves are in calibrated i?-band magnitudes. Epochs in parentheses 
were rejected as outliers (see Section l3.2.1ll . Only a portion of this table is shown here 
to demonstrate its form and content. 


Table 2 

SDSS J0924+0219 Ultraviolet (F275W) Light Curves 


HJD - 2,450,000 


5253.8 

5279.3 

5311.1 

5338.9 

5364.4 

5476.3 
5626.8 

5659.3 

5704.4 

5727.1 


QSO A 


18.730 ±0.001 
19.530 ±0.013 
19.346 ±0.026 
19.184 ±0.055 
19.134 ±0.001 
20.115 ±0.001 
19.104 ±0.030 
19.103 ±0.008 
19.073 ±0.001 
19.108 ±0.001 


QSO B 


20.385 ± 0.016 
21.088 ±0.001 
20.849 ± 0.028 
20.921 ±0.001 
20.706 ±0.015 
21.561 ±0.060 
20.770 ±0.001 
20.644 ± 0.001 
20.526 ± 0.004 
20.644 ± 0.009 


QSO C 


21.612 ±0.001 
22.386 ± 0.005 
22.028 ±0.025 
21.998 ±0.040 
21.952 ±0.016 
22.895 ± 0.023 
21.911 ±0.023 
22.028 ±0.012 
21.812 ±0.029 
21.866 ±0.028 


QSO D 


22.045 ± 0.066 
22.760 ± 0.008 
22.529 ±0.001 
22.498 ± 0.002 
22.360 ±0.032 
23.260 ±0.007 
22.382 ±0.002 
22.397 ±0.026 
22.279 ±0.001 
22.393 ±0.048 


Note.— Light curves are in AB magnitudes. In this filter, the offset from ST to AB 
magnitudes is rriAB — msT = 1.532. 
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Fig. 2.— i?-band light curves for images A+D summed (top curve), B (middle curve), and C (bottom curve) in SDSS 
J0924+0219. There are 191 epochs of SMARTS data, and 60 epochs of Euler data (beginning on day HJD — 2450000 = 5510 
in 2010 and ending on HJD - 2450000 = 6364). 


optical/UV flux ratios. We divided the full ob¬ 
served frame X-ray spectral range (0.4-8.0 keV) 
into soft (0.4-f.3 keV) and hard (f.3-8.0 keV) en¬ 
ergy bands. This corresponds to a full-band range 
of 1-20 keV in the quasar rest frame. We analyzed 
the soft and hard X-ray variability separately to 
test for consistency and to explore any possible 
energy structure in the X-ray continuum source. 

Figure [4] shows the difference UV and X-ray 
light curves between (merged) images A-l-D and 
B, A-l-D and C, and between B and C in the 
top, middle, and bottom panels, respectively. The 
amplitude of the optical microlensing is clearly 
largest during the first season. During the sev¬ 
enth and eighth seasons where we have multi¬ 
wavelength data, the uncorrelated X-ray variabil¬ 
ity is pronounced with respect to the optical and 
UV, while the uncorrelated UV variability is pro¬ 
nounced with respect to the optical. 


3. Analysis 


We analyze the optical, UV, and X - ray li ght 
curves using the technique of Kochane^ (12004) as 
adapted for multiwavelength analyses ( Mor£an_et_aJj 


2008 , 2012 ; Dai et a0201oi: Poindexter fc Kochanek 


20101 ). For the present analysis, we are uninter¬ 


ested in the surface brightness profile of the X-ray 
source, but rather the physical extent, so we ex¬ 
ploit the fact that microlensing variability can 
constrain the half-light radius independent of the 
surface brightness profil e adopted for the sourc e 
(Mortonson et al.l 


s profil e adopted tor the sourc e 
] I 2 OO 5 I : Congdon et al. 20071) . 


Since SDSS J0924-I-0219 has negligible time delays 
(< 10 days) between the im ages as compared t o 
the microlensing timescales ( Morgan et al l liooel) . 
we attribute any residual variability in the dif¬ 
ference light curve between two images to mi¬ 
crolensing rather than intrinsic quasar variability 
(unless the timescale for intrinsic variability is 
shorter than either the time delay or the sam¬ 
pling timescale, which may be important in the 
X-ray data). Our analysis began with creating a 
sequence of strong lens models to characterize the 
macroscopic properties of the lens galaxy. Then, 
using these lens models, we generated detailed 
microlensing magnification patterns to represent 
the distribution of stars in the lens galaxy near 
each of the quasar images, as described in Sec¬ 
tion o We employed a Monte Carlo simulation 
to generate many realizations of the model mi¬ 
crolensing light curves over a range in source size, 
using physically-motivated values of the tangential 
velocity. Finally, we performed a (joint) Bayesian 
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Fig. 3.— Top-left panel: HST UV light curves in AB magnitudes for the f our qu asar images in SDSS J0924-I-0219. Top-right 
panel: full-band X-ray light curves for SDSS J0924-I-0219 from IChen et aO ll2012ri . The remaining panels show the X-ray data 
split into hard and soft bands for each image individually. 


analysis on the acceptable model light curves us¬ 
ing the optical, UV, and X-ray monitoring data 
(Sections ElUHSXl- 

3.1. Lens Models 

We started by modeling the strong lensing in 
SDSS J0924-I-0219 to obtain convergence k and 
shear 7 estimates in the vicinity of each quasar 


image. We modeled the l ens galaxy with LENS- 
MODE jd^ ( Keeton _ 2001 b selecting concentric 
NEW ( Navarro et al.lll996[) and de Vaucouleurs 
mass profiles to represent the dark matter and 
stellar mass components, respectively. Our con¬ 
straints on the lens astrometry are based on the 
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Fig. 4.— Difference light curves in the i?-band (gray dots), UV (blue asterisks), and full X-ray band (black diamonds) for 
SDSS J0924+0219. Note that while the X-ray and UV fluxes for images A and D are summed here, they are treated as separate 
in the analysis. 


HST i?-band i mage from CAS TLES, the details 
of which are in Morgan et al.l ( 20061) . We mod¬ 
eled the le ns as an early-type gal axy at redshift 
zi = 0.39 ( Eigenbrod et al. 20(161) with effective 
radius Ve = 0'.'31±0"02, axis ratio q = 0.92 ±0.02, 
and major axis position angle 9e = —27°±8° (east 
of north). The Einstein ring curve of the quasar 
host galaxy was included as a constraint on the 
mass models ( Morgan et al ]l2006l: l iKochanek et al 


200 ll) . and an external shear was included to ac¬ 
count for any perturbations from other objects 
near the lens or along the line of sight. We 
did not attempt to use the observed flux ra¬ 
tios as constraints on the strong lensing model, 
since they are anomalous due to microlensing 


and millilensing from substructure (IKeeton et al 


2006 : Eigenbrod et al. 200^ Morgan et al. 2006t 


Faure et al.ll201ll : iBate et al.ll201ll) . 

We parametrize the unknown stellar mass frac¬ 
tion k^/k by creating a sequence of ten models 
in which the mass of the de Vaucouleurs compo¬ 
nent is set to the fraction /m/l of the mass of the 
constant mass-to-light model with no dark halo. 
We used a uniform sampling from /m/l = 1-0 and 
0.1 in steps of 0.1, reoptimizing the lens model for 
each case. The convergence and shear estimates 
are provided in Tabled Finally, by adopting a sin¬ 
gular isothermal sphere profile for the lens galaxy, 
we estimated the stellar velocity dispersion in the 
lens to be 215 km s“^. Such estimates are known 
to agree well with dynamical measurements (e.g., 
Bolton et al. 2008h . 


















































































Table 3 

Macroscopic Lens Mass Models 


Jm/l 

A 

Convergence k, 
B C 

D 

A 

Shear 7 

B C 

D 

A 

k*Ik, 

B C 

D 

0.1 

0.73 

0.72 

0.77 

0.74 

0.24 

0.21 

0.29 

0.29 

0.019 

0.017 

0.023 

0.020 

0.2 

0.67 

0.65 

0.70 

0.67 

0.30 

0.26 

0.37 

0.37 

0.041 

0.038 

0.048 

0.041 

0.3 

0.59 

0.58 

0.64 

0.60 

0.36 

0.31 

0.45 

0.44 

0.065 

0.060 

0.083 

0.069 

0.4 

0.53 

0.52 

0.57 

0.53 

0.42 

0.36 

0.53 

0.52 

0.10 

0.09 

0.12 

0.10 

0.5 

0.46 

0.45 

0.51 

0.47 

0.48 

0.40 

0.62 

0.60 

0.14 

0.13 

0.17 

0.14 

0.6 

0.39 

0.38 

0.44 

0.40 

0.54 

0.45 

0.70 

0.68 

0.20 

0.19 

0.24 

0.20 

0.7 

0.32 

0.31 

0.37 

0.33 

0.60 

0.50 

0.78 

0.75 

0.27 

0.26 

0.33 

0.28 

0.8 

0.26 

0.25 

0.31 

0.27 

0.65 

0.54 

0.86 

0.83 

0.38 

0.37 

0.45 

0.40 

0.9 

0.19 

0.18 

0.24 

0.20 

0.71 

0.59 

0.95 

0.90 

0.58 

0.57 

0.65 

0.60 

1.0 

0.12 

0.11 

0.18 

0.13 

0.77 

0.64 

1.03 

0.98 

1.00 

1.00 

1.00 

1.00 


Note.—C onvergence re, shear 7 and the fraction of the total surface density composed of stars re*/re 
at each image location for the series of macroscopic mass models. 


3.2. Monte Carlo Microlensing Analysis 

For each of the 10 macromodels parametrized 
by /m/L: ■we generated 40 unique sets of periodic 
microlensing magnification patterns fo r each im¬ 
age us ing the P^M method described in Kochanel] 
( 20041) ■ The magnification patterns are 8192 x 
8192 pixel arrays with an outer scale radius of 
20i?_B, where Re = 5.7 x {M /cm is 

the Einstein radius of a star of mass M in the 
source plane, yielding a source plane pixel scale of 
1.4 X 1O14(M/M0)1/2 cm. For a black hole with 
mass Mbh = 2.8 x IO^Mq (as estimated from the 
Mgll emission line width in Morgan et al. 2006), 
this pixel scale is 3.4 times the estimated size of 
the black hole, Vg = GMbh/c^- These patterns 


have t wice the resolution of those in I Morgan et al.l 
(20061 ) . which is important for simultaneously re¬ 
solving the optical and X-ray sources. 

Our simulations are carried out in Einstein 
units, where source sizes fg and velocities Ve are 
related to their physical scales by the mass of 
an average lens galaxy star (r^ = 

Ve = Ve{M To convert from Einstein 
units to physical units, we convolved the proba¬ 
bility density for the scaled source size P{fs) with 
the probability density for the scaled effective ve¬ 
locity, P{ve), and a statistical model (i.e., a prior) 
for the true effective source velocity, P{ve)- We 
constructed Pjv^ ) using the method described in 
Kochanei] ( 20041) . adopting the peculiar velocity 


estimates for the redshifts of SDSS J0924-I-0219 
an d the lens galaxy from the m odels presented 
in Mosauera fc Kochanek ( 201ll) . Our velocity 
model includes the 176 km s“^ projected veloc¬ 
ity of the CMB dipole onto the plane of the lens, 
a probability distribution for the one-dimensional 
peculiar velocity dispersion of galaxies at zi with 
an rms of 272 km s“^, and a one-dimensional 
stellar velocity dispersion in the lens galaxy of 
215 km s“^, based on the lens models. An al¬ 
ternate method of converting to physical units is 
to assume a prior on (M). In Section 13. 2. 11 we de¬ 
scribe the optical-only analysis using such a prior. 
The joint optical-LTV analysis is presented in Sec¬ 
tion 13.2.21 and the joint optical-X-ray analysis is 
presented in Section [3.2.31 


3.2.1. Optical Microlensing Models 

The first step in our dual Monte Carlo anal¬ 
ysis was to generate 4 x 10"^ trial light curves for 
each magnification pattern. We performed 40 sep¬ 
arate Monte Carlo simulations of each magnifica¬ 
tion pattern, attempting a grand total of 1.6 x 10^ 
trials for each source size and image. Trials were 
made for a grid of 15 different source sizes span¬ 
ning 13. 5 < log (fs,opt/cm) < 17, yielding 2.4 x 10® 
total trial fits to the i?-band data. We added a 
systematic error of 0.01 mag in quadrature to the 
photometric uncertainties for images A-l-D and B, 
and a systematic error of 0.025 mag for C. This ad- 
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justment yields a reasonable number of fits with 
X^/Ndof ~ 1 and helps compensate for any prob¬ 
lems in the photometry such as flux sharing be¬ 
tween images due to poor seeing. Since we did not 
use the observed flux ratios to constrain our strong 
lens models, we allowed for 0.5 mag of systematic 
uncertainty in the absolute flux ratios relative to 
the macro model. This helps compensate for any 
problems in the macro model magnifications and 
any effects of substructure. The final goodness 
of fit, in the sense that the probability of the data 
given the model is P{D\M) oc e~^ is evaluated 
from the of the light curves, where all data have 
their true statistical weights. While a model mi- 
crolensing light curve is generated for each of the 
four images, the summed model light curve for im¬ 
ages A and D is compared to the summed data for 
these images when computing the x^- We dropped 
all fits with x^/Ndof > 3.0, as they contribute neg¬ 
ligibly to the Bayesian integrals, leaving a total of 
2.8 X 10^ complete fits to the optical light curves. 

In the process of fitting the microlensing vari¬ 
ability, we also identified 18 outliers in the optical 
data due to poor photometry. We considered five 
of the best microlensing fits to the original data in 
which we have first subtracted our model for the 
intrinsic variability. For each example fit, we iden¬ 
tified all observations which were > Scr discrepant 
from adjacent epochs for at least two of the three 
quasar images. These observations had systemat¬ 
ically worse seeing than average. Eighteen obser¬ 
vations were flagged as outliers in all five exam¬ 
ple fits, and we exclude these points from the final 
analysis. These data are enclosed in parentheses in 
Table m It is important to exclude these outliers 
from the analysis because they would otherwise 
dominate the x^ of the fits and significantly reduce 
the number of actual good fits that are retained 
in the Monte Carlo analysis. Given the individ¬ 
ual trials, we then marginalize over the nuisance 
parameters using uniform priors for the trajectory 
starting points, the velocity prior described above, 
and log arithmic priors on the sizes (for full de¬ 
tails. see lKochanekl[2003 IPoindexter fc Kochanek 
2 ninli . 

Figure [5] shows example fits to the observed op¬ 
tical flux ratios. We express the optical source 
radius as the scale radius r^^opt at which the disk 
temperature corresponds to the rest-frame wave¬ 
length of the R band (where T = hv/k). The cen- 
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Fig. 5.— Example optical microlensing model light 
curves (solid curves) for images A-l-D, B, and C. The data 
points with error bars show the data after subtracting the 
intrinsic source variability. The x^ /Ndof is 1-1 for both 
panels. 


tral hole in the disk emissivity is neglected to avoid 
the introduction of an additional parameter, and is 
small enough to be unimportant. From the optical 
data alone, we find a most probable source scale 
radius of Vg^opt = 2 x 10^® cm with a Icr confidence 
interval of 7 x 10^'^ cm < rg ,opt < 4 X 10^*^ cm. 
If we impose a uniform prior on the mean mi¬ 
crolens mass of 0.1 < (M/Mq) < 1.0, we ob¬ 
tain rs,opt = 0.9l();g X 10^® cm. Assuming a typ¬ 
ical inclination of 60° and converting to a half- 
light radius (riy 2 = 2.44rs), we obtain riy 2 ,opt = 


- 1 - 7.7 

- 4.5 


X 10^° cm. When using a uniform prior 


6.9 

on th e mean microlens mass as in [Morgan et al 
( 2006l l. we obtain ri/ 2 ,opt = 3.2 
consistent with their estimates. 


) 4 - 4.8 
^- 2.3 


X 10 
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3.2.2. Optical-UV Joint Analysis 


A 


Following IPai et aL ( 20ir)[) . all physical vari¬ 
ables from the good optical fits (e.g., trajectory 
and effective transverse velocity) were saved for 
the combined optical and UV analysis, and fits to 
the UV light curves were then attempted using 
the same solution but for 12 different UV source 
sizes spanning 13 < log (fg^uv/cm) < 15.75. For 
the UV data, we added a systematic error of 0.01 
mag in quadrature to the photometric uncertain¬ 
ties for each image. This is probably an underes¬ 
timate of the systematic uncertainties in the UV 
fluxes (e.g., as created by intrinsic variability), but 
it allowed us to find reasonable numbers of solu¬ 
tions. When assessing the quality of our fits to 
the observed UV flux ratios, we again allowed for 
0.5 mag of systematic uncertainty in the flux ra¬ 
tios of the macro models. All fits to the UV data 
with x^/Ndof < 30 were retained (note: the joint 
x'^/^dof remains < 3). Figure [5] shows an example 
fit to the observed UV flux ratios. 


The analysis of the combined UV and optical 
solutions yielded a set of joint probability densi¬ 
ties for the variables of interest, namely the UV 
and optical source sizes. The resulting constraints 
on rs_uv and fs,opt are listed in Table IH and the 
probability distributions are shown in Figure [T] 
In physical units, after converting to a half-light 
radius and correcting both sizes for inclination i 
assuming (cosi) = 1/2, the resulting UV size is 
10^“^ cm < ri/ 2 ,uv < 3 x 10^^ cm, while the optical 
size is larger at 3 x 10^® cm < ri/ 2 ,opt < 10^® cm. 
Here, we are adopting the results from using a stel¬ 
lar mass prior, as the probability distributions are 
more stable when altering the limit on the UV 
model light curves. The probability distribution 
for the optical-UV size ratio is shown in Figure [8] 
The probability distribution for /m/l resulting 
from the joint analysis is flat, indicating no prefer¬ 
ence for the stellar mass fraction, and the best-fit 
mean stellar mass is log {{M/Mq)) = —0.5 ± 1.0. 


3.2.3. Optical-X-ray Joint Analysis 

Similar to the optical-UV analysis, all physi¬ 
cal variables from the good optical fits were saved 
for the combined optical and X-ray analysis. The 
same source sizes used in the UV analysis were 
attempted for the X-ray fits, and we analyze the 
microlensing variability in the full, soft, and hard 
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Fig. 6.— Example microlensing model light curves (solid 
curves) for the UV data (top panel), soft X-rays (middle 
panel), and hard X-rays (bottom panel). The data points 
with error bars show the observed data after subtracting 
our model for the intrinsic source variability. 
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X-ray bands separately. For the X-ray data, we 
allowed a systematic uncertainty of 0.1 mag in 
the photometry for all images, as this yielded a 
sufficient number of fits, and we again allowed a 
systematic uncertainty of 0.5 mag in the flux ra¬ 
tios. The former uncertainty should help reduce 
any contamination from intrinsic variability, which 
likely has a shorter timescale than the time delays 
(6 hrs in the observed frame for X-ray emission 
within a radius of ^ 7rg in SDSS J0924-I-0219). 
All fits with x^/A^dof < 2.4, 1.5, and 1.76 for the 
full, soft, and hard bands, respectively, were re¬ 
tained. Figure [6] shows example X-ray fits from 
our joint analysis in the soft and hard bands. 

The hard-band size falls right at the pixel scale 
of our magnification patterns, so we can only re¬ 
port an upper limit in this case, while the soft- 
band size is resolved (Figure [7]). Since the sizes 
resulting from the joint analysis of the R and soft 
bands are the most tightly constrained, we quote 
the best optical size from this analysis. Hereafter, 
we report the results from using a stellar mass 
prior since they are more stable against the cho¬ 
sen cutoff in x^, as with the UV results. Table |4] 
lists the resulting sizes. 


In physical units, after converting to a half- 
light radius and assuming a mean inclination 
of (cosi) = 1/2, the resulting optical size is 
larger than the p reviously reported result 
in Morgan et al.l (l2006ll . The hard-band analysis 


results in a la upper limit of ri/ 2 .hard < 10^® cm 
(or ri/2,hard < 6 X 10^^ Cm at 95% confidence). 
Compared to the optical size, the best-fit soft- 
band size (3 X 10^^ cm) is 13 times smaller. The 
probability distributions for the ratios of the soft- 
and hard-band sizes to the optical size are shown 
in Figure 13 and the individual probability densi¬ 
ties are shown in Figure |3 


We find a best-fit mean stellar mass of 
log ((M/Mq)) = — 0.9to'9: where (M) is calcu¬ 
lated by combining the scaled effective source ve¬ 
locity distribution with the prior probability dis¬ 
tribution for the true source velocity. The width 
of the (M) distribution shown in Figure [T3 re¬ 
flects the uncertainty in the effective velocity since 
{M) (X Vg'^. Finally, the probability distribution 
for /m/l resulting from our joint analysis is flat, 
indicating no preferenc e for the stellar mass frac¬ 
tion, in agreement with Morgan et al.l ( 2006ll . 



(<M/Mg>'/^ cm) 



(<M/Mg>'/^ cm) 


Fig. 7.— Probability distributions for the scaled UV 
^s,uv X-ray (red points), and the optical rs,opt 
(black points), from the joint analysis (top: optical-UV, 
middle: optical—soft-X-ray, bottom: optical-hard-X-ray). 
The distributions have been scaled to a unit peak height 
to ease visual comparisons. The dashed line shows the es¬ 
timate including a prior of 0.1 < (M/Mq) < 1.0 on the 
mass of the stars. The vertical lines show the Einstein ra¬ 
dius of the average-mass star and the pixel scale of the 
magnification map in the source plane. The expected Vg is 
at 4.1 X (M/Mq >1/2 cm. 
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aAs.R 

Fig. 8.— Probability distributions for the ratios of 
source sizes (dotted: optical—UV, solid: optical-soft X-ray, 
dashed: optical-hard X-ray), scaled to a unit peak height 
to ease visual comparisons. No sizes have been corrected for 
inclination. For a given ratio, the two sizes do not necessar¬ 
ily result from the same trial trajectory. Therefore, these 
distributions do not take into account the covariances be¬ 
tween the source sizes and may be broader than the joint 
probability distributions. 


4. Discussion 


is significantly larger, by l.Scr, than the theo¬ 
retical thin-disk estimate based on the observed, 
magnification-corrected HSTI-han d flux, sim i¬ 


lar to the results for other systems (iPoolev et al 


2007 Morganet^^ 2010: Blackburn e et ahl 2011 


Hainline et alJ l2013t iMosauera et alJ 120131) . The 
accretion disk radius based on our optical-only 
analysis is con sistent with bu t slight ly larger than 


the result in Morgan et alJ (120061 ). while that 


based on the joint analysis is shifted to even larger 
sizes. Therefore, the combination of using merged 
A and D images in the i?-band and including mul¬ 
tiwavelength data has likely affected the resulting 
optical size. Since images A and D could not 
be reliably separated in the ground-based data, 
systematic errors m ay have biased the result in 
Morgan et al.l ( 2006ll . who trie d to treat the image 


fluxes separately. Note that Flovd et al. ( 2009h 


found an upper bound on the optic al radius that 
was m uch larger than the result in Morgan et al.l 
(|2006l) when adopting the same priors on (M). 
Aside from merging the images, based on our 
analysis of the optical, UV, and X-ray emitting 
regions, the inclusion of UV and X-ray monitor¬ 
ing data seems to provide significant additional 
information, leading to a more robust constraint 
on the optical size and further challenging the 
standard thin disk model. 


We have improved upon and extended the mi- 
crolensing analysis of SDSS J0924-I-0219 by in¬ 
cluding more recent i?-band monitoring and more 
stringent constraints on the microlensing models 
provided by simultaneous space-based UV and X- 
ray monitoring. After attempting over 100 million 
Monte Carlo fits to the light curves, we marginal¬ 
ize over the nuisance variables (trajectory starting 
points and directions) to generate joint probabil¬ 
ity densities for the variables of interest {ve , (M ), 

K*, Ts). 

Adopting the results from the joint analy¬ 
sis, we find a Icr confidence interval for the 
inclination-corrected optical half-light radius of 
1.8 X lO^^cm < ri/2,opt < f.4 x 10^®, approxi¬ 
mately Icr larger than the previous measurement 

of 6.4 X 10^"^c m _ < ri/2,opt < 2.3 x 10^® cm 

( Morgan et al.l l2006l ). When scaled to a rest 
wavelength of 2500A and converted to the scale 
radius where the photon energy matches the 
disk temperature (T = hu/k)^ this size corre¬ 
sponds to logi? 25 oo/cm = 15.2to;3. This size 


SDSS J0924-I-0219 exhibits pronounced mi¬ 
crolensing variability in the UV compared to the 
optical, indicating a comparatively small extent 
for the UV continuum emitting region, as expected 
for an accretion disk. We find a Icr confidence 
interval of 10^^ cm < ri/ 2 ,uv < 3 x 10^^ cm, 
which is significantly smaller (by a factor 8) than 
the optical radius. If we assume that the UV 
flux originates from the thermal portion of the 
disk, our joint optical-UV analysis yields a tem¬ 
perature profile slope of /3 = 0.46 ± 0.46, where 
T oc r~^. While the most probable /3 estimate 
from our analysis suggests a shallower temper¬ 
ature profile than expected for a standard thin 
disk, we cannot rule out the standard thin disk 
expectation of /3 = 0.75 given our uncertainties. 
Howeve r, our result is mildly incons istent with 


that of I.Timenez-Vicente et al.l (j20141 ). who find 


an average slope of ,5 = 1.33 ± 0.36 in a sample 
of eight leri s ed qu asars. For SDSS J0924-I-0219, 
Flovd et ‘ail ( 2009ll obtain j3 = 1.dtp's based on 


the anomalous A/D flux ratio in the optical and 
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Table 4 

Resulting Sizes for SDSS J0924+0219 


Analysis 

Band 

log fs [(M/M©)^/^ cm] 

logri /2 [cm] 

Opt.- 

-Only 


15.2t°j 


Opt. 

-UV 

Opt. 


15.8t°i 



UV 

14.5t?:^ 

14.9t°:® 

Opt. 

-Soft 

Opt. 





Soft 

14.3t°j 


Opt.- 

-Hard 

Opt. 





Hard 

< 14.8 (68%) 

< 15.1 (68%) 

Opt. 

-Full 

Opt. 

1 q+0.3 

-LO.O_o 4 

15.4toi 



Full 


14.8t°j 


Note.— All size estimates result from imposing a prior on the mean 
stellar mass of 0.1 < (M/Mq) < 1.0. The optical and UV rj ^/2 estimates 
have been corrected assuming (cosi) = 1/2. 


near-IR bands at a single epoch. However, this 
estimate depends critically on the absence of sub¬ 
structure and on the choice of logarithmic or lin¬ 
ear priors on the source size. Therefore, we can- 
not make a clear c omparison, as the analysis of 
Keeton et al. ( 2006l l strongly suggests that at least 


part of the anomalous A/D flux ratio is not due 
to microlensing. 

The microlensing variability in SDSS J0924-I-0219 
is even more pronounced in the X-rays compared 
to the optical and UV, indicating an even smaller 
extent for the X-ray continuum emitting region. 
From the soft-band joint analysis, the most proba¬ 
ble size ratio is Copt/rsoft = 13. We were unable to 
constrain the lower limit for the hard-band source 
size due to the sparse sampling of the X-ray data 
and the pixel scale of the magnification patterns, 
but we find an upper limit of r’i/ 2 ,hard < 10^^ cm 
at 68% confidence. These results confirm that 
the X-ray continuum emitting regions in quasars 
are compact, supporting models involving an X- 
ray emitting region concentrated near the inner 


for six other svstems (e.g.. Morgan et al. 

2008 

2012: Chartas et al.l 

2 OI 2 I: Mosauera et al. 

2013 

Blackburne et al. 2013. 20141). Furthermore, while 


the three X-ray size measurements are formally 
consistent, our results suggest that the soft X- 


ray-emitting region is more extended than the 
hard X-ray region, although still very compact 
with respect to the accretion disk. 

Note that the optical and UV half-light radii 
have been corrected assuming an average disk in¬ 
clination of (cos/) = 1/2, and the sizes can be 
rescaled as ri /2 oc (cos/)“^/^. However, other 
sources of error in physical parameters such as 
the black hole mass and projected velocity should 
not have an appreciable effect on our estimated 
source sizes. The uncertainty in the black hole 
mass would only affect the size of the central hole 
in the disk. This would have little effect in the op¬ 
tical and the UV, where the disk model is applica¬ 
ble, and where the effective size is relatively large 
compared to the last circular orbit. Uncertainties 
in the projected velocities are fully included and 
have less effect on the results than naive expec¬ 
tations, as discussed in Kochanek (2004). Lastly, 
any additional unvarying flux at the location of 
each image, such as that due to the Einstein ring 
visible in the HST H-band image, will not affect 
our results, as they use only the time-varying flux 
ratios between the quasar images. 

Assuming that microlen sing is the sol e cause 
of the flux ratio anomaly, iMorgan et al.l ( 2006h 
predicted that the image A/D magnitude differ¬ 
ence {rriA — ttid) would increase by at least 0.25 
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Fig. 9.— Probability densities for the half-light radii of the optical, UV, and X-ray emitting regions, scaled to a unit peak 
height for ease of visual comparison. In the top panel, the soft X-ray, hard X-ray, and UV probability densities are shown as a 
solid blue curve, a dashed red curve, and a dotted black curve, respectively. In the bottom panel, the distributions from our 
optical-only microlensing analysis are shown in red and are broader than those resulting from our joint (soft-band) analysis, 
shown in blue. The dashed curves in the bottom panel show the result of imposing a uniform prior on the mean microlens 
mass of 0.1 < M/Mq < 1.0. All optical and UV sizes have been corrected for a disk inclination of 60°. The vertical line 
shows the Schwarzschild radius Rbh = ‘IGMbhI of a 2.8 x 10 ®Mq black hole dMorgan et al.|[2006ll . The last stable orbit 
for a Schwarzschild black hole is at SRbh- The small artifacts in the curves result from coarse binning in the integration of 
probabilities. 


mag over 3 years with a 70% likelihood. How¬ 
ever, the flux of image D has remained on average 
at least a factor 12.5 fainter than that of image 
A in the optical, and on average a factor 20 ± 1 
(9±3) fainter in the UV (full-band X-rays). There¬ 
fore, we still cannot rule out a significant influence 
from substructure in the lensing potential. How¬ 
ever, since we allow a 0.5 mag uncertainty in the 
(un-microlensed) flux ratios of the quasar images, 
our results should not be sensitive to the effects of 
millilensing. Furthermore, by using a time series 
analysis, our microlensing models do not rely on 


precise flux ratio measurements but rather on the 
time signature of m icrolensing (IKochanekl 12004 ; 


Morgan et al. 2012ll . Future high-resolution mid- 


IR observations of this system, such as from the 
James Webb Space Telescope, should help con¬ 
strain the effect of substructure in the lens galaxy 
by providing the intrinsic flux ratios. 

The simulations presented here could be im¬ 
proved by using dynamic magnification pat- 


terns (i.e., recalculated 


Poindexter fc Kochanek 


for e ach 


(120131) , and iBlackburne et al 


epoch) 


as m 


Mosauera et al.l 

( 20141) ■ although 
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Fig. 10.— Probability distribution for the mean mi¬ 
crolens mass (M) for the lens galaxy resulting from our 
joint optical—soft-X-ray analysis. 


the results should not change significantly. Our 
results would also be improved by running a joint 
optical-UV-X-ray analysis, but that is too com¬ 
putationally expensive for our present facilities. 
Additionally, note that our model for the intrinsic 
variability is simply based on the correlated vari¬ 
ability between the quasar images, and our anal¬ 
ysis could be improved by comparing this model 
to what we expect based o n observed light curve s 
for non-lensed quasars (e.g.. lMacLeod et^l2010ll . 
Finally, we have neglected the gravitational leas¬ 
ing effects of the black hole which may be im¬ 
portant on such small scales, so our simulations 
may be improv ed by adopting a fully relativistic 
disk model (see Chen et al JEll). We leave these 
improvements for a future study. 

With further X-ray monitoring for more sys¬ 
tems, future studies will be able to more precisely 
explore the energy structure of the quasar X-ray 
continuum emission region and test correlations 
such as the one between black hole mass and 
X-ray size (e.g., iMosouera et al.M2013ll . Contin¬ 
ued analysis of the hard and soft X-ray bands 
should improve our constraints and help distin¬ 
guish betwee n physical model s , such as the X- 
ray model of Gardner fc Done ( 2014ll . and those 


involving the presence of reflection components 
(Fabian et al. 2004; see review by Uttley et al. 
2014). For example, a more extended soft X- 
ray compone nt would be suppo r ted b y the physi¬ 
cal models of Gardner Sz Done ( 2014ll . who show 
that propagation-plus-refiection models fail to ex¬ 
plain the soft X-ray lag behavior in the AGN PG 
1244-1-026. As magneto-hydrodynamical accretion 
disk models invol ving an X-ray corona continue 


to improve (e.g., iNoble et al. 


I 2 OO 7 : Penna et al 


2010; iSchnittman et al.l l2006t iSchnittman et al 
2013 ; IPexter fc BlaesI 2014ll . so will our under¬ 


standing of the nature of the various X-ray struc¬ 
tures observed through microlensing. With the 
advent of new surv ey data such as Pan-STARRS 
( Kaiser^^_ah 2010) an d the Dark Energy Survey 
(|Honscheid et al.l 120081) . new lenses will become 
available for a more extensive analysis of quasar 
continuum e mitting regions. T he Large Synoptic 
Sky Survey ( Ivezic et al. 20081) is expected to un¬ 
cover thousands of new lensed quasars, providing 
well-sampled light curves with multiwavelength 
microlensing information, which will improve our 
understanding of the temperature profile of ac¬ 
cretion disks and the source of various types of 
emission in AGN. 


This material is based upon work supported 
by the National Science Foundation (NSF) under 
grant AST-1211146 to CWM. We also gratefully 
acknowledge support by Ghandra X-Ray Center 
award 1170050. CSK is supported by NSF grant 
AST-1009756, and FC and GM are supported by 
the Swiss National Science Foundation (SNSF). 
This research made extensive use of the USNA 
high performance computing cluster. We thank 
the reviewer for valuable comments which im¬ 
proved the manuscript. 

Facilities: CTIO:2MASS (ANDICAM), HST 
(NICMOS, ACS), Eulerl.2m, CXO (ACIS) 

REFERENCES 

Bate, N. F., Floyd, D. J. E., Webster, R. L., & 
Wyithe, J. S. B. 2008, MNRAS, 391, 1955 

—. 2011, ApJ, 731, 71 

Blackburne, J. A., Kochanek, C. S., Chen, B., 
Dai, X., & Chartas, G. 2013, ArXiv e-prints, 
arXiv: 1304.1620 


16 




























































—. 2014, ApJ, 789, 125 

Blackburne, J. A., Pooley, D., & Rappaport, S. 

2006, ApJ, 640, 569 

Blackburne, J. A., Pooley, D., Rappaport, S., & 
Schechter, P. L. 2011, ApJ, 729, 34 

Bolton, A. S., Treu, T., Koopmans, L. V. E., et al. 
2008, ApJ, 684, 248 

Cackett, E. M., Fabian, A. C., Zogbhi, A., et al. 

2013, ApJ, 764, L9 

Chang, K., & Refsdal, S. 1979, Nature, 282, 561 

Chartas, G., Agol, E., Eracleous, M., et al. 2002, 
ApJ, 568, 509 

Chartas, G., Kochanek, C. S., Dai, X., et al. 2012, 
ApJ, 757, 137 

Chen, B., Dai, X., Baron, E., & Kantowski, R. 

2013, ApJ, 769, 131 

Chen, B., Dai, X., Kochanek, C. S., et al. 2011, 
ApJ, 740, L34 

—. 2012, ApJ, 755, 24 

Congdon, A. B., Keeton, C. R., & Osmer, S. J. 

2007, MNRAS, 376, 263 

Corrigan, R. T., Irwin, M. J., Arnaud, J., et al. 
1991, AJ, 102, 34 

Dai, X., Chartas, G., Agol, E., Bautz, M. W., & 
Garmire, G. P. 2003, ApJ, 589, 100 

Dai, X., Kochanek, C. S., Chartas, G., et al. 2010, 
ApJ, 709, 278 

DePoy, D. L., Atwood, B., Belville, S. R., et al. 
2003, in Society of Photo-Optical Instrumenta¬ 
tion Engineers (SPIE) Conference Series, Vol. 
4841, Instrument Design and Performance for 
Optical/Infrared Ground-based Telescopes, ed. 
M. lye & A. F. M. Moorwood, 827-838 

Dexter, J., & Blaes, O. 2014, MNRAS, 438, 3352 

Eigenbrod, A., Courbin, F., Dye, S., et al. 2006, 
A&A, 451, 747 

Eigenbrod, A., Courbin, F., Meylan, G., et al. 
2008a, A&A, 490, 933 


Eigenbrod, A., Courbin, F., Sluse, D., Meylan, G., 
& Agol, E. 2008b, A&A, 480, 647 

Eulaers, E., Tewes, M., Magain, P., et al. 2013, 
A&A, 553, A121 

Fabian, A. C., Rees, M. J., Stella, L., & White, 
N. E. 1989, MNRAS, 238, 729 

Fame, C., Sluse, D., Cantale, N., et al. 2011, A&A, 
536, A29 

Floyd, D. J. E., Bate, N. F., & Webster, R. L. 

2009, MNRAS, 398, 233 

Gardner, E., & Done, C. 2014, MNRAS, 442, 2456 

Garmire, G. P., Bautz, M. W., Ford, P. G., 
Nousek, J. A., & Ricker, Jr., G. R. 2003, in 
Society of Photo-Optical Instrumentation En¬ 
gineers (SPIE) Conference Series, Vol. 4851, 
Society of Photo-Optical Instrumentation En¬ 
gineers (SPIE) Conference Series, ed. J. E. 
Truemper & H. D. Tananbaum, 28-44 

Goicoechea, L. J., Shalyapin, V. N., Koptelova, 
E., et al. 2008, New A, 13, 182 

Hainline, L. J., Morgan, C. W., MacLeod, C. L., 
et al. 2013, ApJ, 774, 69 

Honscheid, K., DePoy, D. L., & for the 

DES Collaboration. 2008, ArXiv e-prints, 
arXiv:0810.3600 

Inada, N., Becker, R. H., Buries, S., et al. 2003, 
AJ, 126, 666 

Ivezic, Z., Tyson, J. A., Acosta, E., et al. 2008, 
ArXiv e-prints, arXiv:0805.2366 

Jimenez-Vicente, J., Mediavilla, E., Kochanek, 
C. S., et al. 2014, ApJ, 783, 47 

Kaiser, N., Burgett, W., Chambers, K., et al. 

2010, in Society of Photo-Optical Instrumen¬ 
tation Engineers (SPIE) Conference Series, Vol. 
7733, Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series 

Keeton, C. R. 2001, ArXiv Astrophysics e-prints, 
arXiv:astro-ph/0102340 

Keeton, C. R., Buries, S., Schechter, P. L., & 
Wambsganss, J. 2006, ApJ, 639, 1 


17 



Keeton, C. R., Gaudi, B. S., & Fetters, A. O. 2005, 
ApJ, 635, 35 

Kochanek, C. S. 2004, ApJ, 605, 58 

Kochanek, C. S., & Dalai, N. 2004, ApJ, 610, 69 

Kochanek, C. S., Keeton, C. R., & McLeod, B. A. 
2001, ApJ, 547, 50 

Kochanek, C. S., Morgan, N. D., Falco, E. E., 
et al. 2006, ApJ, 640, 47 

Krist, J. E., Hook, R. N., & Stoehr, F. 2011, in 
Society of Photo-Optical Instrumentation En¬ 
gineers (SPIE) Conference Series, Vol. 8127, 
Society of Photo-Optical Instrumentation En¬ 
gineers (SPIE) Conference Series 

Laor, A. 1991, ApJ, 376, 90 

Lehar, J., Falco, E. E., Kochanek, C. S., et al. 
2000, ApJ, 536, 584 

Lewis, G. F., & Irwin, M. J. 1995, MNRAS, 276, 
103 

MacLeod, C. L., Ivezic, Z., Kochanek, C. S., et al. 

2010, ApJ, 721, 1014 

Magain, P., Courbin, F., & Sohy, S. 1998, ApJ, 
494, 472 

McLeod, B. A., Bernstein, G. M., Rieke, M. J., & 
Weedman, D. W. 1998, AJ, 115, 1377 

Mediavilla, E., Munoz, J. A., Kochanek, C. S., 
et al. 2011, ApJ, 730, 16 

Morgan, C. W., Kochanek, C. S., Dai, X., Morgan, 
N. D., & Falco, E. E. 2008, ApJ, 689, 755 

Morgan, C. W., Kochanek, C. S., Morgan, N. D., 
& Falco, E. E. 2006, ApJ, 647, 874 

—. 2010, ApJ, 712, 1129 

Morgan, C. W., Hainline, L. J., Chen, B., et al. 
2012, ApJ, 756, 52 

Morgan, N. D., Chartas, G., Malm, M., et al. 2001, 
ApJ, 555, 1 

Mortonson, M. J., Schechter, P. L., & Wambs- 
ganss, J. 2005, ApJ, 628, 594 

Mosquera, A. M., & Kochanek, C. S. 2011, ApJ, 
738, 96 


Mosquera, A. M., Kochanek, C. S., Chen, B., et al. 
2013, ApJ, 769, 53 

Mosquera, A. M., Munoz, J. A., Mediavilla, E., & 
Kochanek, C. S. 2011, ApJ, 728, 145 

Munoz, J. A., Mediavilla, E., Kochanek, C. S., 
Falco, E. E., & Mosquera, A. M. 2011, ApJ, 
742, 67 

Navarro, J. F., Frenk, C. S., & White, S. D. M. 
1996, ApJ, 462, 563 

Noble, S. C., Leung, P. K., Gammie, C. F., & 
Book, L. G. 2007, Classical and Quantum Grav¬ 
ity, 24, 259 

Peeples, M. S., Schechter, P. L., & Wambsganss, 
J. K. 2004, in Bulletin of the American Astro¬ 
nomical Society, Vol. 36, American Astronomi¬ 
cal Society Meeting Abstracts, 1392 

Penna, R. F., McKinney, J. C., Narayan, R., et al. 
2010, MNRAS, 408, 752 

Peterson, B. M., Alloin, D., Axon, D., et al. 1992, 
ApJ, 392, 470 

Poindexter, S., & Kochanek, C. S. 2010, ApJ, 712, 
658 

Poindexter, S., Morgan, N., & Kochanek, C. S. 

2008, ApJ, 673, 34 

Pooley, D., Blackburne, J. A., Rappaport, S., & 
Schechter, P. L. 2007, ApJ, 661, 19 

Pooley, D., Blackburne, J. A., Rappaport, S., 
Schechter, P. L., & Fong, W.-f. 2006, ApJ, 648, 
67 

Pooley, D., Rappaport, S., Blackburne, J. A., 
Schechter, P. L., & Wambsganss, J. 2012, ApJ, 
744, 111 

Rathna Kumar, S., Tewes, M., Stalin, C. S., et al. 
2013, A&A, 557, A44 

Schnittman, J. D., Krolik, J. H., & Hawley, J. F. 
2006, ApJ, 651, 1031 

Schnittman, J. D., Krolik, J. H., & Noble, S. C. 
2013, The Astrophysical Journal, 769, 156 

Shakura, N. L, & Sunyaev, R. A. 1973, A&A, 24, 
337 


18 



Shappee, B. J., Prieto, J. L., Grupe, D., et al. 
2014, ApJ, 788, 48 

Smith, H. J., & Hoffleit, D. 1963, Nature, 198, 650 

Tewes, M., Courbin, F., Meylan, G., et al. 2013, 
A&A, 556, A22 

Teyssier, R., Moore, B., Martizzi, D., Dubois, Y., 
& Mayer, L. 2011, MNRAS, 414, 195 

Vaughan, S., Edelson, R., Warwick, R. S., & Utt- 
ley, P. 2003, MNRAS, 345, 1271 

Wambsganss, J., Witt, H. J., & Schneider, P. 1992, 
A&A, 258, 591 


This 2-column preprint was prepared with the AAS lAT^ 
macros v5.2. 


19 



